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— ^The  results  calculated  predict  a significant  change  of  the  wave  pattern 
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The  present  method  can  be  applied  for  prediction  of  the  bow  wave.  The 
numerical  calculation  is  performed  for  a double-wedge-shaped  model  used 
by  Standing  (1974).  The  calculated  bow  wave  pattern  shows  good  agreement 
with  the  corresponding  experimental  result. 
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ABSTRACT 


The  problem  of  a potential  flow  around  a surface- 
piercing strut  at  a high  Froude  number  is  proposed.  The 
strut  can  be  of  asymmetric  cross-section  and  is  not  neces- 
sarily slender.  The  exact  boundary-value  problem  is  formu- 
lated and  is  approximated  by  a mathematical  model  which 
satisfies  the  linearized  free-surface  boundary  conditions 
on  the  undisturbed  free-surface  and  the  rigid-body  boundary 
condition  on  the  actual  surface  of  the  strut. 

It  is  found  that  the  problem  can  be  broken  into  two 
parts,  one  of  which  is  a conventional  infinite-fluid  problem, 
the  other  a problem  which  can  be  formulated  as  a two-dimen- 
sional unsteady-flow  problem. 

A numerical  scheme  is  developed  and  applied  for  three 
cases,  namely,  aspect  ratios  of  0.5,  1.0  and  3.0  . The 
cross-section  of  the  strut  is  a lens  shape,  the  thickness/ 
length  ratio  of  which  is  0.1  . Side-force  coefficients  and 
wave  patterns  are  calculated  for  various  Froude  number. 

The  side-force  coefficient  calculated  for  0.5  aspect 
ratio  shows  good  agreement  with  the  experimental  result  for 
a flat  plate  of  the  same  aspect  ratio  and  angle  of  attack. 

The  side-force  coefficient  for  aspect  ratios  1.0  and  3.0  is 
shown  to  change  with  respect  to  Froude  number  in  a pattern 
similar  to  that  for  aspect  ratio  0.5  . Also,  it  is  shown 
to  approach  the  analytical  limit  values  as  the  Froude  number 
increases. 
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The  results  calculated  predict  a significant  change  of 
the  wave  pattern  on  the  suction  side  of  the  strut  as  the 
aspect  ratio  changes. 

The  present  method  can  be  applied  for  prediction  of  the 
bow  waves.  The  numerical  calculation  is  performed  for  a 
double-wedge-shaped  model  used  by  Standing  (1974).  The 
calculated  bow  wave  pattern  shows  good  agreement  with  the 
corresponding  experimental  result. 
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The  potential  flow  around  a surface-piercing  strut 
travelling  at  a high  Froude  number  is  the  subject  of  inves- 
tigation in  this  work.  The  cross  section  of  the  strut  can 
be  of  either  symmetric  or  asymmetric  form.  The  solutions 
are  limited  to  thin  struts  with  a small  degree  of  lateral 
asymmetry.  Both  the  leading  edge  and  trailing  edge  should 
be  sharp.  These  restrictions  are  necessary  so  that  the 
waves  generated  by  the  presence  of  the  strut  are  small  and 
so  the  free-surface  boundary  condition  may  be  linearized. 

It  was  probably  in  a case  of  surface-piercing  hydro- 
foils (Nishiyama,  1964)  that  the  problem  involving  lateral 
asymmetry  was  studied  seriously  for  the  first  time.  How- 
ever attempts  in  this  direction  did  not  go  beyond  lifting- 
line theory.  Newman  comprehensively  discussed  the  problem 
of  asymmetrical  surface-piercing  struts  after  the  fashion 
of  lifting-surface  theory.  Kern  (1973)  also  used  lifting- 
surface  theory  and  solved  the  inverse  problem.  Daoud 

i 

• \ 

(1973)  succeeded  to  solve  the  proposed  problem  for  the  case  i 

of  a yawed  surface-piercing  flat  plate.  However  Daoud 's  \ 

work  shows  how  complicated  it  is  to  solve  the  two-dimen-  | 

sional  singular  integral  equation  of  linear  lifting-surface 
theory  with  the  presence  of  the  free  surface,  for  even  the 
simplest  case.  i 

In  the  present  work,  the  boundary-value  problem  is 
solved  following  a method  which  is  related  to  slender-ship 

theory.  j 

{ 

t 
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Tuck  (1963)  introduced  the  method  of  matched  asymptotic 
expansion  into  ship  hydrodynamics  to  develop  his  slender- 
ship  theory.  Tuck's  method  made  it  possible  to  take  full 
advantage  of  the  slender-ship  characteristics.  That  is, 
the  two-dimensional  approximation  and  the  rigid-wall  free- 
surface  condition  can  be  applied  near  the  ship  provided  the 
Froude  number  is  order  of  unity.  However,  the  wave  resis- 
tance according  to  Tuck's  theory  becomes  infinite  except 
for  a ship  with  a cusped  wedge  bow  or  a pointed  bow. 

Ogilvie  (1972)  assumed  the  longitudinal  perturbation 
near  the  bow  of  a slender  ship  (bow  near  field)  to  be  still 
small  compared  with  the  transverse  perturbation  but  much 
larger  than  it  is  near  the  center  part  of  the  ship.  With 
this  assumption  he  obtained  a bow  near  field  where  the 
linearized  free-surface  condition  applied  instead  of  the 
rigid-wall  condition. 

The  Laplace  equation  remained  two-dimensional.  He 
also  used  an  assumption  that  the  solution  is  independent  of 
downstream  conditions.  With  these  assumptions,  he  solved 
the  problem  in  the  bow  near  field.  Later,  Ogilvie  (1975) 
succeeded  in  obtaining  the  far-field  solution  which  matched 
with  the  bow-near-field  solution  near  the  ends  and  matched 
with  the  usual  near-field  solution  near  the  middle  part  of 
the  body.  When  he  applied  this  modification,  the  wave 
resistance  from  slender-ship  theory  exactly  coincided  with 
the  Michel-integral  value.  His  assumption  that  the  longi- 
tudinal perturbation  is  much  larger  in  the  bow  near  field 


Miiiiiitfililliaiiai 


3 


than  it  is  near  the  center  part  of  the  ship  can  be  under- 
stood in  that  the  Froude  number  in  the  bow  near  field  is 
much  higher  than  the  usual  Froude  number,  which  depends 
on  ship  length.  This  is  because  the  Froude  number  can  be 
taken  as  a measure  of  the  relative  magnitude  of  inertial 
forces  with  respect  to  gravitational  forces  in  the  interior 
of  a fluid  region.  His  other  assumption,  stating  that  the 
solution  is  independent  of  downstream  conditions,  can  be 
understood  more  clearly  by  changing  the  longitudinal  axis 
into  time  domain;  X = UT.  Then  the  problem  is  reduced  to 
a two-dimensional  unsteady  flow  problem  with  the  presence 
of  the  free  surface.  Hence  Ogilvie's  method,  which  depends 
on  those  two  assumptions,  may  be  considered  equivalent  to 
a time  stepping  method.  Ogilvie  solved  the  problem  with  an 
approximation  of  a thin  body  or  thin  bow.  This  is  an  extra 
simplification  to  the  problem  which  is  not  necessary  in 
the  slender-ship  theory.  Daoud  (1975)  solved  the  same 
problem  for  the  case  of  a bow  with  finite  thickness.  He 
satisfied  the  body  boundary  condition  on  the  actual  surface 
of  the  ship.  Hence  his  method  includes  the  diffraction  of 
diverging  waves  by  the  body.  The  time-stepping  method  can 
be  applied  to  problems  involving  lateral  asymmetry.  Chapman 
(1976)  solved  the  problem  for  the  case  of  a yawed  slender 
surface-piercing  flat  plate. 

However,  since  the  time-stepping  method  reduces  the 
problem  to  a two-dimensional  unsteady  flow  problem,  an 
initial  condition  at  the  bow  is  necessary.  Ogilvie,  Daoud, 
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and  Chapman  all  assumed  that  there  is  no  disturbance  ahead 
of  the  bow  or  leading  edge.  This  condition  may  be  reason- 
able for  a really  slender  body.  However  even  for  fairly 
slender  bodies  such  as  fine-ship  models,  it  is  observed 
that  the  wave  elevation  at  the  bow  is  of  an  appreciable 
magnitude  (see  Standing  (1974)  for  example).  Also,  the 
acceleration  becomes  infinite  at  the  leading  edge  according 
to  the  assumption  of  no  disturbance  ahead  of  the  leading 
edge.  In  the  bow  near  region,  the  inertia  force  is  expected 
to  dominate  the  gravitational  force,  but  it  is  not  infinite. 
Because  of  these  reasons  the  effects  of  the  disturbance 
ahead  of  the  leading  edge  are  expected  to  be  significant 
even  for  a fairly  slender  body  and  by  no  means  should  be 
neglected  in  the  strut  problem.  Furthermore,  Ogilvie's 
method  depends  on  the  assumption  of  the  slenderness  of  the 
body  and  so  cannot  be  applied  to  a strut  which  is  not 
slender . 

Chapman  showed  that  the  nondimensional  problem  of  the 
flow  around  a yawed  surface-piercing  flat  plate  can  be 
determined  by  only  two  nondimensional  parameters  according 
to  the  time  stepping  method.  One  of  them,  named  n t is 
the  product  of  the  Froude  number  and  the  square  root  of  the 
aspect  ratio,  while  the  other,  called  pn  , is  the  product 
of  the  Froude  number  and  the  ratio  of  angle  of  attack  to 
the  square  root  of  the  aspect  ratio.  On  the  other  hand, 
it  is  known  that  the  free-surface  condition  can  be  approxi- 
mated by  the  pressure-relief  surface  condition  at  a very 
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high  Froude  number.  Because  of  this  reason  it  can  be 
concluded  that  the  flow  around  a yawed  surface-piercing 
flat  plate  should  be  independent  of  Froude  number  and  so 
independent  of  n or  yg  at  a high  Froude  number.  Thus 
there  are  significant  upstream  effects  of  any  disturbance, 
and  so  it  seems  that  the  assumption  of  the  slenderness  of 
the  body,  which  leads  to  the  time  stepping  method,  is  too 
strong.  There  is  an  essential  difficulty  in  including 
effects  of  aspect  ratio,  at  least  at  a high  Froude  number. 

Ogilvie  (1975)  proposed  an  idea  to  overcome  this 
difficulty.  If  the  problem  could  be  broken  into  two  parts, 
one  of  which  is  a conventional  infinite-fluid  problem,  the 
other  a problem  in  wave  propagation  to  which  can  be  reason- 
ably applied  the  time  stepping  method,  the  problem  can  be 
solved  without  great  difficulty.  The  proposed  method  can 
be  referred  to  as  a hybrid  method,  since  the  problem  will 
be  solved  partially  three-dimensionally  and  partially  two- 
dimensionally . 

This  idea  of  the  hybrid  method  is  developed  and  applied 
to  the  problem  of  the  surface-piercing  strut  at  a high 
Froude  number.  The  side  force  acting  on  the  strut  and  the 
wave  pattern  around  the  strut  are  calculated  numerically 
for  different  Froude  numbers  and  aspect  ratios. 
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II.  PROBLEM  FORMULATION  AND  METHOD  OF  SOLUTION 
Formulation  of  the  problem 

The  strut  is  assumed  to  be  fixed  in  an  incident  stream 
of  velocity  U.  The  coordinate  system,  as  shown  in  Figure  1, 
has  the  2-axis  in  the  upward  direction,  the  Y-axis  positive 
to  starboard  and  the  X-axis  parallel  to  and  in  the  same 
direction  as  the  incident  stream.  The  origin  of  the  coor- 
dinate system  is  set  at  the  intersection  of  the  trailing 
edge  of  the  strut  and  the  undisturbed  free  surface. 


Figure  1.  The  coordinate  system 


The  chord  length  and  the  depth  of  the  strut  are  repre- 
sented by  L and  D . Also  H,  SS,  SP,  SB,  LE  and  TE  will 
be  used  ns  abbreviated  notations  of  the  whole  surface,  the 
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starboard  side  surface,  the  port  side  surface,  the  bottom 
surface,  the  leading  edge  and  the  trailing  edge  of  the 
strut  respectively.  The  equation  of  the  hull  is  given  as: 

H(X,Y,Z)  = 0 . 

Notation  "H"  will  be  used  to  denote  either  the  surface  of 
the  strut  or  the  function  defining  the  surface  of  the  strut. 

The  free  surface  is  given  as: 

Z = ^d(X,Y)  . 

The  strut  is  assumed  to  be  thin,  with  a small  degree 
of  lateral  asymmetry.  Also,  the  shape  and  dimensions  of 
cross  sections  must  vary  slowly  in  the  X direction. 

In  case  of  an  asymmetric  strut,  a wake  exists  behind 
the  strut.  We  assume  that  the  wake  occupies  a thin  sheet 
that  joins  the  strut  at  its  trailing  edge  and  extends 
infinitely  downstream,  the  surface  of  the  sheet  being  given 
as 

Y = W(X,Z)  . 

The  usual  assumptions  are  made  about  the  fluid  being 
inviscid  and  incompressible.  The  fluid  motion  is  assumed 
to  be  irrotational  except  on  the  wake.  There  exists  a 
velocity  potential,  which  we  write  in  the  following  form: 

UX  + 4>(X,Y,Z)  . 

This  velocity  potential  satisfies  the  Laplace  equation  in 
the  fluid  domain  except  at  the  wake: 

^ ^ ^ ° ' 
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the  dynamic  free-surface  boundary  condition  on  the  free 
surface : 

g^D  + + ^2  2U<^}  * 0 on  Z = Cd(X,Y)  , 

the  kinematic  free-surface  boundary  condition: 

(U+$x)^Dx  ^Y^Dy  - <I'z  = 0 on  Z = Cd(X,Y), 

and  the  rigid-surface  boundary  condition  on  the  surface  of 
the  strut: 

= -Unx  on  H, 

where  n(X,Y,Z)  is  a unit  normal  vector,  positive  into  the 
fluid  domain.  nx  , ny  and  n2  are  the  components  of  n in 
X,  Y,  and  Z directions.  The  dynamic  condition  on  the  wake, 
which  states  that  there  is  no  discontinuity  of  the  pressure 
across  the  wake,  is, 

(2U$x+  = (2U<I>x  + $4  + + ^l) 

^ ^ ^ Y=W+0  A X Y Z 

(1) 

and  the  kinematic  condition  which  states  that  wake  consists 
of  stream  lines  is, 

4>Y  - {{U+$x)Wx  + ^'zWzJ  = 0 on  Y = W . (2) 

Finally,  we  must  have  a radiation  condition  which  states 
that  the  wave  propagation  dies  off  faster  upstream  than 
downstream  of  the  strut. 

The  method  of  solution 


^ The  exact  solution  of  the  boundary- value 
problem  as  specified  cannot  be  determined  by  existing 


analytic  techniques.  Since  our  case  includes  asymmetric 
struts,  wo  have  to  face  the  infinite  perturbation  velocity 
along  the  sharp  leading  edge  as  well  as  extra  conditions  at 
the  wake  in  addition  to  the  usual  difficulties  common  to 
the  problem  of  a symmetric  body  on  the  free  surface.  In 
particular,  the  location  of  the  wake  is  not  given  explicitly. 
Also  a discontinuity  of  the  wave  elevation  is  observed  at 
the  trailing  edge.  The  precise  investigation  of  the  last 
problem,  i.e.,  the  discontinuity  of  the  wave  elevation  at 
the  trailing  edge  and  the  phenomena  associated  with  it  will 
be  left  for  future  works.  However  our  theory  should  include 
the  discontinuity  of  the  wave  level  at  the  trailing  edge 
and  should  be  able  to  predict  the  fluid  motion  downstream 
of  the  trailing  edge  in  some  consistent  way.  This  will  be 
discussed  later.  Similarly  the  infinite  perturbation 
velocity  at  the  leading  edge  should  be  included  in  our 
theory. 

The  major  difficulties  left  are  the  nonlinearity  of 
the  conditions  at  the  free  surface  and  the  wake.  As  a first 
step  toward  making  it  possible  to  solve  the  problem,  we 
shall  approximate  these  nonlinear  conditions  with  linear 
conditions.  The  free-surface  conditions  may  be  linearized 
by  expanding  the  potential  in  a Taylor  series  about  Z=0 


and  discarding  all  high  order  terms.  Thus  we  obtain: 
g^D  + U$x  =0  on  Z=0 

and 


on 


Z=0 
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or,  combining  them  together, 

<l>xx  ® Z=0, 

where 

K = g/U^  . 

The  theoretical  determination  of  the  exact  position  of  the 
wake  is  beyond  our  resources  even  in  an  infinite  fluid.  In 
this  work,  the  wake  is  assumed  to  be  a plane  sheet  parallel 
to  the  direction  of  the  incident  flow.  We  define  as 

the  plane  sheet  that  the  wake  occupies.  Then  is 

given  as: 

Y=0  , X>0  , 0>Z>-D. 

We  also  define  as  the  starboard  side  surface  of 

and  Sy,”  as  the  port  side  surface  of  Sy,  . These  defini- 
tions are  shown  in  Figure  2.  This  model  of  the  wake 
violates  the  kinematic  condition  at  the  wake.  The  kinematic 
condition  at  the  wake,  (2) , should  be  replaced  by  a weaker 
condition  which  states  the  continuity  of  the  velocity  in 
the  normal  direction  across  the  wake.  Let  ny,  be  the 
normal  vector  on  Sy,^  the  surface  of  the  wake,  and  A^^iy, 
be  the  difference  of  the  velocity  in  the  normal  direction 
across  the  wake.  Then  the  condition  stated  above  can  be 
given  as 

^^nw  ~ ‘^nw  4.  “ "^nw  _ = 0 . 

Sw  Syi 


The  equation  (1),  the  dynamic  condition  at  the  wake,  becomes 
after  linearization: 


Figure  2.  The  wake  model 


Hence  the  potential  jump  across  the  wake 

= {$  - $ } 

Sw 


is  independent  of  X.  Let  'I's(Q)  be  the  value  of  $ on 
the  surface  of  the  strut.  Also  define 
as  follows: 

lim  (Q)  “ ^st  '■ 

Q-+TE  on  SS 

and 

lim  “ ’^’pt^^)  • 

Q->TE  on  SP 

Using  these  notations,  the  condition  obtained  above  can  be 
stated  as 


Z 
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The  value 

4>st(Z)  - 'J>pt(Z) 

is  the  circulation  around  the  strut,  by  definition;  hence, 
finally, 

A$(Z)  = $st(Z)  - <I>pt(Z)  = rD(Z)  on  S„  . 

For  more  discussion  concerning  the  wake  and  the  conditions 
at  the  wake  we  refer  to  Robinson  & Laurmann  (1956,  1.15). 

The  Laplace  equation,  the  body  boundary  condition  and 
the  radiation  condition  remain  unchanged. 

It  is  convenient  to  rewrite  the  defining  equations  in 
nondimensional  form.  Nondimensional  coordinates  are: 

X = X/L,  y = Y/L,  and  z = Z/L  , 
and  the  nondimensional  velocity  potential  is 
((I  = ^/UL  ; 

also,  the  nondimensional  wave  height  is 

^ ' 

Other  nondimensional  values  are  as  following: 

AR(=  aspect  ratio)  = D/L 
Fn(=  Froude  number)  = U//gL 
r(z)  = rj3(Z)/UL  . 

The  unit  normal  vector,  n(x,y,z)  , its  components,  njj,ny, 
n2  and  the  wake,  Sy,  , and  its  surfaces,  Sy,'*'  and  Sy,"  , 
are  to  be  defined  in  nondimensional  space.  Then  the 
nondimensional  boundary-value  problem  is  given  as  follows: 
The  Laplace  equation  in  the  fluid  domain: 

't’xx  "^yy  "^zz  “ 


0, 


(3) 
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the  dynamic  free  surface  condition: 
r,  + = 0 on  z=0  , 

and  the  kinematic  free-surface  condition; 


Cv  = 


on  z=0 


(or  the  combined  linearized  free-surface  condition: 
4>xx  ^n^'Pz  ~ ® z=0)  , 

the  body  boundary  condition: 

(|)n  = -njj  on  H, 

the  condition  on  the  wake: 

^’i’nw  ~ ‘J’nwl 


■'nw 


=’w 


= 0 


on  S„  , 


A(J)(2)  — (f)g^(z)  — (^p^(z)  — r(z)  on 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


step  2 Although  the  boundary  value  problem  we  have 
to  solve  is  linear  now,  it  is  still  an  extremely  difficult 
one  to  solve.  Therefore  we  will  follow  the  idea  of  the 
hybrid  method  to  further  facilitate  finding  a solution  of 
the  problem. 

Now  we  break  the  velocity  potential  <)>  into  two  parts 
as  follows: 

= 4>o  + <<>i  • 

Since  the  ({(g  portion  of  the  problem  is  to  be  solved 
three-dimensional ly,  the  free-surface  condition  should  be 
simple.  There  are  two  choices.  One  of  them  is  the 
condition 


4*0  = 0 on  z=0  , 


and  the  other  is  the  condition 

(|)(,^  = 0 on  z=0  . 

To  decide  which  condition  to  choose,  we  refer  to  the 
beautiful  pictures  issued  by  NSRDC  (Rothblum,  Mayer  and 
Daily,  1972)  of  yawed  sharp-edged  surface-piercing  struts 
moving  at  high  speed.  In  many  of  those  pictures,  cavitation 
is  observed  along  the  leading  edge  on  the  suction  side. 

This  coincides  with  the  assumed  infinite  perturbation 
velocity  along  the  leading  edge.  However  it  is  clearly 
shown  in  those  pictures  that  the  cavitation  region  is 
clearly  separated  from  the  free  surface.  Also  the  free- 
surface  elevation  near  the  leading  edge  is  quite  different 
from  what  would  be  obtained  by  assuming  that  the  flow  on 
the  free  surface  rounds  the  sharp  leading  edge.  Therefore 
we  propose  a condition  stating  that  the  flow  on  the  free 
surface  does  not  round  the  sharp  leading  edge.  To  be 
consistent  with  this  condition,  we  must  use  the  former 
free-surface  condition, 

(})|j  = 0 on  z=0  , (10) 

instead  of  the  rigid-wall  free-surface  condition. 

On  the  surface  of  the  strut  we  use  the  condition: 

<t)n  = "Hv  on  H,  (11) 

“n 

for  the  (t)^  portion  of  the  velocity  potential.  Other 
conditions  which  should  be  satisfied  by  the  ({ij  portion 
of  the  velocity  potential  are  the  Laplace  equation  in  the 
fluid  domain  (except  at  the  walce)  , 


4>o  + <J>0  ‘J’o  = 0 » 

®xx  yy  “zz 


the  conditions  at  the  wake. 


A<t)„  = 0 


on  Sw  f 


A4)o  (z)  = r,  (z) 
Also,  a condition. 


on  St 


lim  (()q  = 0 

X-*— <» 


is  assumed,  so  that  A can  be  determined  uniquely. 

0 

Because  the  Aj  portion  of  the  velocity  potential  is 
defined  by 

<|)  1 = - 4>o  ' 

it  should  satisfy  the  Laplace  equation  in  the  fluid 
domain  (except  at  the  wake) : 


A + Ai  + 4>i  = 0 , 

^xx  yy  zz 


the  linearized  free-surface  condition 
j.  I I*— = n <-> 


<^1  + Fn  <t>i  + Fn^Ao  =0  on  z=0 

XX  z z 


(or,  separated  into  the  dynamic  and  the  kinematic  free- 
surface  conditions. 


= 0 

on 

z=0 , 

r « = (t)  + 

^ “z 

on 

o 

II 

N 

the 

condition  on 

the 

surface 

of  the 

o 

II 

c 

■e- 

on 

H, 

and 

the  conditions  at 

the  wake: 

AA , =0 

nw 

on 

Swf 
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A4)i  (z)  = r{z)  - r^Cz)  = Tj  (z) 


on  S„. 


A radiation  condition  is  also  necessary.  However,  since 
the  time  stepping  method  will  be  used  eventually,  and  hence 
the  waves  will  propagate  downstream  only,  it  is  not 
necessary  to  define  the  radiation  condition  in  terms  of  (Pg 
and  <()j  . The  condition  for  the  uniqueness  of  the  solution 
will  suffice  for  a condition  at  infinity. 

It  is  observed  that  at  a high  Froude  number  the  first 
crest  of  the  transverse  waves  lies  far  behind  the  strut 
and  the  diverging  waves  are  dominant  around  and  close 
behind  the  strut.  The  crests  of  these  diverging  waves  are 
almost  parallel  to  the  incident  stream  (or  the  x-axis) ; 
hence  it  is  reasonable  to  assume  that 


=o  U-  A 
9x  l3y  ' dz 


for  the  (|)j  portion  of  the  problem.  This  is  because  these 
waves  motions  are  represented  by  the  ({)  j portion  of  the 
velocity  potential.  Hence  the  (t>j  portion  of  the  problem 
can  be  reasonably  approximated  by  a sequence  of  two-dimen- 
sional problems  in  the  transverse  planes,  the  potentials 
for  which  satisfy  the  two-dimensional  Laplace  equation  in 
the  fluid  domain,  except  at  , 

({i  +(tij=0,  (20) 

‘yy  zz 

the  linearized  free  surface  condition 

+ ^^n^<l>n  “0  on  z=0  (21) 

‘xx  z “z 


(or,  separated, 
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r.  + Fp^(|),^  = 0 on  z = 0,  (22) 

and 

Cx  * <{>  1 + <{>0  on  z=0)  , (23) 

z z 

the  conditions  at  the  wake; 

A({),  = 0 on  Sw  , (24) 

^nw 

A(J)i(x)  = rj(z)  on  S„  . (25) 

The  body  boundary  condition  on  the  surface  of  the  strut  is 
to  be  satisfied  approximately  in  the  two  dimensions  of  the 
transverse  planes  as: 

4)i^=0  onH/  (26) 

where  N is  the  two-dimensional  normal.  A condition 
stating  that 

())j  ->  0 as  X 

or  /y 2+z 2 -y  oo  (27) 

is  assumed  in  order  to  make  the  solution  unique. 

Equation  (17)  , the  free  surface  condition  for  ())j  , 
shows  that  (J)j  is  the  velocity  potential  for  a fluid 
motion  which  could  be  caused  by  an  applied  pressure  distri- 
bution on  the  free  surface.  If  we  let  the  equivalent 
pressure  distribution,  nondimensionalized  by  pU^  , be 
represented  by  PQ(x,y)  , its  form  is  given  by: 

P0x(5t,y)  = Fn^<})02  ‘ 

The  disturbance  caused  by  a concentrated  pressure  distri- 
bution on  a stream  is  known  to  die  off  much  more  quickly 
upstream  than  downstream  (see,  for  example,  Peters  1949). 
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Because  of  this  reason,  it  seems  reasonable  to  assume  that 
the  <j)j  portion  of  the  problem  is  independent  of  the 
condition  downstream.  This  condition  can  be  understood 
more  clearly  y transferring  the  x-axis  into  a nondimen- 
sional  time  domain.  Then  the  (})j  portion  of  the  problem 
can  be  formulated  as  a two-dimensional  unsteady-flow 
problem.  Especially  behind  the  trailing  edge  where  there 
is  no  rigid  boundary,  the  4)1  problem  can  be  considered  as 
a generalized  two-dimensional  Cauchy-Poisson  problem. 

Hence  the  solution  of  the  4)i  problem  in  the  region  behind 
the  trailing  edge  can  be  obtained  just  by  evaluating  a 
double  integral  (for  example,  see  Lamb,  1932) . The 
notation  x and  t will  be  freely  used  interchangeably 
for  convenience  hereafter. 

The  (J)p  portion  of  the  problem  corresponds  to  the 
case  of  infinite  Froude  number  and  so  includes  no  effect 
of  gravity.  In  the  following  chapters,  the  <}>o  portion  of 
the  problem  will  be  mentioned  as  the  "zero-gravity  problem" 
or  simply  as  the  "4)g  problem".  The  <))j  portion  of  the 
problem  includes  all  the  wave  motion  and  hence  will  be 
mentioned  as  the  "wave  problem"  or  as  the  "(})j  problem". 
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III.  SOLUTION  OF  THE  ZERO-GRAVITY  PROBLEM  (({)„  PROBLEM) 

The  problem  of  the  (J)^  portion  of  the  velocity  poten- 
tial can  be  looked  at  as  a three-dimensional  problem  of  a 
body  in  an  infinite  fluid.  Methods  originally  developed 
for  a lifting  or  nonlifting  body  problem  in  an  infinite 
fluid  (Hess  & Smith,  1966;  Hess,  1972;  Webster,  1975;  etc.) 
can  be  modified  so  that  they  can  be  applied  to  the 
problem. 

However,  a rather  basic  method  is  to  be  used  in  this 
work.  Define  a control  surface  S^  and  a coordinate 
system  O^np  as  in  Figure  3,  where 

St  = H + Sf  + Sr  + + S„"  . (28) 


Sy^"*"  and  Sw“  were  previously  defined  as  the  starboard- 
and  port-  side  surface  of  the  wake,  extending  to  infinity 
downstream.  However,  S^'*’  and  Sw~  are  now  redefined  as 
the  parts  of  the  original  Sw"^  or  S^"  bounded  by  Sr 
and  H.  Sr  is  the  part  of  the  hemispherical  surface 
defined  by: 

p < 0 , 

and  bounded  by  Sv,"*"  and  S^~ . H is  the  surface  of  the 
strut  except  for  the  trailing  edge.  Sf  is  the  part  of 
the  undisturbed  free  surface  bounded  by  H,  Sj.  , Sy,"*”  and 
S..”  . Also,  V is  the  fluid  domain  bounded  by  the  closed 
control  surface  S^  • 

According  to  Green's  formula,  4)5  at  a point  P(x,y,z) 


in  V can  be  written  in  terms  of  (J), 


and  ())j  over  S^ 
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ns  Ff)l]ows; 


<J)o(P)  = ;^  ||{<^o  (Q)  •Gn{P;Q)  - 4>o^  (Q) 'G  (P;Q)  }ds  (Q)  , 


where  n is  positive  into  the  fluid.  Choose  G(P;Q)  of 
the  following  form: 

- FTW  ' 

where  P is  the  point  (x,y,z)  and  Q is  (Cfh/P)/  and 
r(P;Q)  = { (x-5) ^+(y-n) 2+(z-p) 

/ (31) 

r'(P;Q)=  { (x-C) ^+(y-n) ^+(z+p) 

Then  G = 0 on  the  undisturbed  free  surface,  and  because 
of  the  free-surface  condition  for  (j);,  , 

(j)^=0  on  z = 0. 

The  contribution  of  the  integral  over  Sf  in  Equation  (29) 
is  zero.  Also  the  contribution  of  the  integral  over  Sj. 
is  known  to  vanish  as  R ->  <»  , from  condition  (15)  , which 

states  that  there  is  no  disturbance  far  upstream  of  the 
strut.  Since  n on  , which  we  shall  denote  by  n^;  , 

equals  -(n)  on  S^”  , the  integral  over  Sy,"*"  and  Sw" 
can  be  combined  as  an  integral  over  Sy,  as  follows: 


r*  + O •" 


1 
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= K ff'**.  «n„  - A+nw  G)  as  . 


Because  of  the  shape  of  and  the  conditions  (13)  and 

(14) , the  above  integral  can  be  reduced  to 


4it 


dpFj  (p)  • Gr,  dC 


Finally,  after  making  R infinite,  equation  (29)  is  reduced 


Gn  - 4.0  (Q)  G}  dS 


■ 

417 


+ ^ dpr„(p)  dC  . 


Because  4.^  over  H is  known  from  the  body  boundary 
condition  and  Fq  given  as  the  potential  jump  at  the 

trailing  edge,  the  only  unknown  in  the  right-hand  side  of 


(32)  is 


on  H , the  surface  of  the  strut.  Now  take  the 


limit  as  P approaches  the  point  Qg  on  the  surface  of 
the  strut: 

lim  r|)^(P)  = j4>o(0s)  + A ff  Uo(Q)Gn  " 4>o  (Q)G}dS 
P->Oc-  ^ J 
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dpi’  (()) 
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The  bar  through  the  integral  signs  denotes  that  a principal 
value  interpretation  should  be  made.  Finally 


^ ]^Uo(Q)Gn  - <Po^(Q)G}  dS 


H 


0 


+ 2?  dpr,(p) 


'n  ; 


(33) 


-AR  0 

this  is  an  integral  equation  to  be  solved. 

Equation  (33)  is  a Fredholm  integral  equation  of  the 
second  kind,  which  can  be  solved  numerically  without 
difficulty.  Then,  using  (to  on  the  surface  of  the  strut 
in  the  equation  (32)  , the  velocity  potential  ()>j  anywhere 
in  V can  be  obtained.  Differentiating  both  sides  of 
(32)  with  respect  to  z,  we  obtain  the  equation: 


<b„  (P)  = 

- 4Tr 


- <J>o„  Gz>  dS 


n 


H 


0 


4Tr 


dpFj  (P) 


Gfiz  * 


(34) 


-AR  0 

Thus  we  can  obtain  the  value  of  on  the  free  surface. 


which  is  necessary  to  solve  (J)^  problem.  Also,  since  the 
(J)j  distribution  on  the  surface  of  the  strut  is  known,  we 
can  calculate  the  pressure  distribution  and  finally  the 
forces  acting  on  the  strut  due  to  the  (J)^  portion  of  the 
problem. 


I 


IV.  THE  WAVE  PROBLEM  ((l)i  problem) 


In  chapter  II,  the  <p  ^ portion  of  the  problem  was 
formulated  as  a two-dimensional  unsteady-flow  problem.  To 
solve  this  problem,  we  use  a method  which  is  essentially 
equivalent  to  what  Chapman  (1976)  used  to  solve  his  problem. 
This  method,  roughly  speaking,  consists  of  two  iterated 
steps.  One  of  the  steps  is  to  obtain  enough  information  to 
solve  the  two-dimensional  boundary-value  problem  at  a 
certain  time  (station)  from  the  solution  at  the  previous 
time  (station) . The  other  one  is  the  solution  of  the  two- 
dimensional  boundary-value  problem  at  each  time  or  station. 

Consider  a problem  at  t = t^  or  at  a corresponding 
transverse  plane  given  as  x = . Let  CF,  CH  and  CR  be 

the  contours  shown  in  Figure  4 and  CT  be  the  closed  control 
contour  given  as: 

CT  = CH  + CF  + CR  ; 

also  let  S be  the  region  bounded  by  CT.  Using  Green's 
formula,  4)  j at  a point  P(y,z)  in  S can  be  written  in 
terms  of  (b , and  d),  on  CT  as  follows: 

4.j(P)  = i {<J),(Q)  Gj,(P;Q)  - (Q)  G(P;Q)}  dZ.  (35) 

We  choose  the  two-dimensional  Green  function,  G(P;Q)  , as: 


G = log 


r^(P;Q)  . 
r(P;Q)  ' 
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Figure  4 . Contours  in  4)  j problem 
where  P is  the  point  (y;z)  and  Q is  (n;P)  and 
r = { (y-n)  ^+(z-p)  r'  = { (y~n)  *+(z+p) 


Then,  for  Q on  the  free  surface, 

G(P;Q)  = 0. 

We  also  have  conditions  for  (|>j: 

4) , = 0 on  CH , 

and 

lim  4>i  0 on  CR  . 

R-^oo 

When  these  conditions  are  applied.  Equation  (35)  is  reduced 
to : 


4>,  (P) 


1 

2ir 


J 4>,(Q) 

CF 


GN(P;Q)d)l 


J 4>i(Q) 

CH 


GN{P;Q)d)l  . 


(36) 


Assume  Qs(>^/y)  is  a point  on  CH  and  take  a limit  as  P 
approaches  Qg  . Then  Equation  (36)  becomes  an  integral 


•■•ff 
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equation. 


(|)i(Qs)  - j ^ I 

CF  CH 


(37) 


which  can  be  solved  for  on  CH  if  the  value  of  (j)  j 

on  the  free  surface,  CF  , is  known.  Hence  the  value  of 
(t)j  on  the  free  surface  is  all  the  information  we  need  for 
solving  the  two-dimensional  problem  at  a certain  time 
(station).  The  value  of  (j)  j on  the  free  surface  can  be 
obtained  from  the  solution  at  the  previous  time  (station) . 
This  can  be  performed  through  the  free-surface  conditions. 
The  linearized  dynamic  free-surface  condition  is 
C + = 0; 

and  the  kinematic  condition  is 
Ct  = <t>iz  + <f>o^  • 

We  may  use  these  conditions  to  obtain  the  value  of 
the  free  surface  at  t = t^  in  the  following  way; 


(p  1 on 


<Pi 


= <pl 

+ Atp  <t>i 

t-tn 

t-tp-Atp 

= <t>i 

- Fp'AtpC 

t-tp-Atp 

t=tn-^tn/2 


t=tp-Atn/2 , 


where ! 


= tp  - tn-l  . 

This  is  because  (|>i^  = • Also  the  value  of  C at 

t-tp-Atj,/2  can  be  obtained  as  follows: 

cl  = C I 

't=tn-Atn/2  ' t=tn-Atn-Atn-i/2 


+ (1/2) (Atn+Atn-i)<  l| 

' t-tn-i 


Since  Ct-  = + <l>n  r we  have  : 

^ 'z  z \ 

{ 

1 

cl  = cl 

't=tn-Atj^/2  ' t=tn-Atn-Atn-i/2 

+ (1/2)  (Atn+Atn-i)  (<}>i  +<}>(,  H • 

2 Z I 4-=:4-  . 

z.  x^n-l 

(38) 

Finally  we  obtain: 

-i’ll  - ^n'^tnfcl 

t=tn  't=tn-Atjj  I-  I t=tn-Atn-Atn-i/2 

+ (1/2)  (Atn+Atn-i)  (<j),  +4>o  )| 

z lt=tn_iJ. 

(39) 

Thus  the  value  of  <()  j on  the  free  surface  and  the  wave 
height  are  obtained  in  terms  of  the  solution  at  the  previous 
time  or  at  the  transverse  section  ahead  of  the  present 
section.  The  accuracy  of  approximation  can  be  made  to  any 
required  degree  by  making  At  smaller.^ 

Now  the  value  of  (fij  on  the  free  surface  is  known. 

Then  the  integral  equation,  (37)  , which  is  a Fredholm 
integral  equation  of  the  second  kind,  can  be  solved  and  the 
value  of  (J)  j on  CH,  the  surface  of  the  strut,  can  be 
obtained.  We  may  differentiate  both  sides  of  Equation 
(36),  to  obtain  the  equation: 

^Chapman  (1976)  showed  that  this  procedure  is  stable 
numerically. 
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G«  (P;Q)dJl  . 


Since  we  have  the  value  of  (J)j  on  CF  and  CH,  the  value  of 

(p,  on  the  free  surface  can  be  calculated  using  Equation 
^ z 

(40) , and  then  it  will  be  used  in  calculating  the  value  of 
(t)j  on  the  free  surface  at  the  next  time  station. 

We  need  an  initial  condition  to  start  the  time-stepping 
procedure  described  above.  The  initial  condition  is  given 
as  the  condition  of  no  disturbance  far  ahead  of  the  leading 
edge  which  can  be  stated  as: 

= 0 and  ^ = 0 at  t = t^  = -T; 

where 

T >>  1. 

This  condition  includes  the  condition: 

d) , =0  on  the  free  surface  at  t = t,  . 


The  time-stepping  procedure  can  be  started  in  the  following 

way  using  these  initial  conditions.  At  t = ti,  (J) , = 0; 

z 

hence 

^t|  " I 

't=ti  ^lt=ti 


Since  C = 0 at  t = t^  , 

rl  = (l/2)At2Ct| 

lt=ti+At2/2  't=ti 

= (l/2)At2<l>o  I 
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Because  of  the  dynamic  free  surface  condition 
<<>1^  = '• 

and  since  <j)  j = 0 at  t = 

(|)i  I = -F-2  (l/2)At2(t)o  I 

It=t2  2lt=t2  • 

Now  that  the  value  of  ({)  j at  t = t2  and  the  value  of  t, 

at  t = t]^+(l/2)At2  are  obtained,  we  can  use  the  usual 

time-stepping  procedure  given  by  the  equations  (38)  and 

(39) , and  we  can  go  on.  It  must  be  noticed  that  in  the 

region  ahead  of  the  leading  edge,  it  is  not  necessary  to 

solve  the  integral  equation.  Since  in  this  region  the 

contour  CH  does  not  exist,  we  can  use  the  equation  (40) 

directly  to  calculate  (})  j on  the  free  surface. 

z 

Just  behind  the  trailing  edge,  the  wave  elevation  is 
not  continuous.  The  two-dimensional  unsteady  fluid  motion 
started  from  a discontinuous  wave  elevation  is  )cnown  to 
be  highly  singular  at  the  point  where  the  discontinuity 
of  the  wave  elevation  was  originally  located  (see  Lamb, 
1932,  for  example).  Therefore  the  solution  is  expected 
to  be  singular  along  the  x-axis  behind  the  trailing  edge. 
However,  for  numerical  solution,  we  can  approximate  the 
discontinuous  wave  elevation  by  a steep  but  continuous 
wave  elevation.  The  time-stepping  procedure  proposed  can 
be  applied  for  this  approximated  wave  elevation  at  the 
trailing  edge  and  the  results  should  be  reasonable  for 
phenomena  the  scale  of  which  is  not  too  small. 


V.  NUMERICAL  RESULTS  AND  CONCLUSIONS 


A computer  program  was  developed  based  on  the  analysis 
given  in  the  previous  chapters  to  solve  the  problem 
numerically.  Figure  5 shows  in  plan  view  the  strut  for 
which  the  numerical  calculation  performed,  and  Figure  6 
shows  the  strut  orientation.  The  results  were  obtained 
for  three  cases  of  aspect  ratios,  namely  0.5,  1 and  3,  at 
different  Froude  numbers.  The  angle  of  attack  was  tan"^ 
0.08  (=5°). 

Some  examples  of  the  results  of  intermediate  steps 
or  the  results  of  (i)^  problem  are  shown  in  Figures  7 - 
Figure  9 . In  Figure  7 is  shown  a typical  spanwise  distri- 
bution of  side  force  in  nondimensional  form  due  to  the 
())g  portion  of  the  potential.  The  side-force  coefficient 
is  plotted  vs.  aspect  ratio  in  Figure  8.  Also  a typical 

example  of  the  distribution  of  4)5  over  the  free  surface 

z 

is  shown  in  Figure  9 . It  should  be  mentioned  that  the 
corresponding  fictitious-pressure  distribution  over  the 
free  surface  which  excites  the  (J)j  portion  of  the  problem 
is  finite  and  smooth  everywhere  on  the  free  surface. 

Final  results  are  shown  in  Figures  10  - 16.  The 
total  side-force  coefficient  for  the  case  of  0.5  aspect 
ratio  plotted  in  Figure  10  shows  good  agreement  with  the 
experimental  result  for  a yawed  surface-piercing  flat 
plate  reported  by  Van  Den  Drug  (1971) , provided  Froude 
number  is  not  too  small. 
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Figure  6 


The  strut  in  the  incident  stream 
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Present  Theory 


Limit  Value  for  Fn=°“ 


Experimental  Data  for  a Rectangular  Flat  Plate 
from  Van  Den  Brug  (1971)  (Plate  D,  AR=0.5) 


Figure  10.  Side-force  coefficient, 
tan  a = 0.08,  AR  = 0.5 
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The  present  theory  does  not  involve  the  assumption  of 
snail  aspect  ratio.  In  Figure  11,  the  total  side-force 
coefficients  are  plotted  vs.  Froude  number  for  different 
aspect  ratios.  For  infinite  Froude  number,  the  linearized 
free-surface  condition  is 

(J)  = 0 on  z = 0 . 

This  is  the  case  of  the  problem;  hence  it  is  reasonable 

that  the  total  side-force  coefficients  are  approaching  the 
values  from  the  4)5  problems  for  corresponding  aspect 
ratio  as  the  Froude  number  increases.  As  the  Froude  number 
decreases  from  infinity,  the  total  side-force  coefficients 
increase  to  the  maximum  value  near  Froude  number  equal  0.4  . 
It  should  be  noticed  that  these  peak  values  are  much 
bigger  than  the  values  from  the  inf inite-Froude-number 
problems.  This  means  that  the  effects  of  the  waves  on  the 
side-force  cannot  be  neglected  unless  the  Froude  number 
is  extremely  high. 

For  Froude-number  equal  zero,  the  linearized  free 
surface  condition  reduces  to 

4)2  = 0 on  z = 0. 

This  is  the  case  of  the  double-body  problem.  The  results 
from  the  double-body  problems  for  corresponding  aspect 
ratios  are  also  shown  in  Figure  11.  However,  the  present 
theory  cannot  be  expected  to  converge  to  these  values  as 
the  Froude  number  becomes  zero,  since  an  explicit  assumption 
of  high  Froude  number  was  made  in  formulating  and  solving 
the  problem.  In  Figures  12  - 16  are  shown  the  wave  patterns 


Figure  13.  Contour  map  of  nondimensional 
wave  height,  AR  = 0.5,  = 1 


igure  16.  Contour  map  of 
nondimensional  wave  heigh 
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created  by  the  strut.  The  diffraction  of  the  waves  by  the 
body  can  be  noticed  in  those  figures.  Another  feature  worth 
noticing  is  that  the  wavf  pattern  changes  much  on  the  suction 
side  of  the  strut  as  the  aspect  ratio  changes  while  it 
changes  little  on  the  pressure  side  of  the  strut.  Also 
the  disturbances  ahead  of  the  leading  edge  increase  with  -I 

i 

the  increase  cf  the  aspect  ratio.  1 

j 

The  present  method  can  be  applied  for  predicting  the  ; 

wave  pattern  near  the  bow  of  a ship. 

Standing  (1974)  measured  wave  patterns  for  a model 

i 

shown  in  Figure  17.  Present  theory  was  applied  also  to  this 


Figure  17 . The  double-wedge  model  in  a plan  view 

case  and  was  compared  with  Standing's  experimental  result 
in  Figure  18.  The  present  theory  shows  good  agreement  with 
experimental  result  in  the  region  very  near  the  bow. 

The  improvement  is  significant  compared  with  the  time- 
stepping method  shown  in  Figure  19,  calculated  by  Daoud 
(1975) . In  the  latter,  disturbances  upstream  (as  given 
in  the  present  ())j  problem)  are  not  considered. 


From  the  results  shown  in  this  chapter,  it  can  be 
concluded  that  the  present  theory  succeeded  in  including 
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Present  theory 


Figure  18.  Contour  map  of  wave  height,  gC/U^,  present 
method  compared  with  measurements  by 
Standing  (1874)  (U=1.83  m/soc,  F,j=0.23) 
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the  disturbance  ahead  of  the  leading  edge  and  the  effects 
of  the  aspect  ratio,  and  hence  it  was  developed  in  the 
right  direction. 


APPENDIX  DOUBLE  BODY  PROBLEM 


The  method  used  for  solving  the  (])„  problem  can  be 
applied  to  a lifting-body  problem  in  an  infinite  fluid. 

It  should  be  noted  that  the  present  method  has  an  advantage 
that  the  solution  of  a lifting-body  problem  which  satisfies 
the  body-boundary  condition  on  the  actual  surface  of  the 
body  can  be  obtained  by  solving  a Fredholm  integral  equation 
of  the  second  kind.  Also  it  can  be  applied  to  a body  with 
sharp  leading  edge.  In  this  appendix,  the  present  method 
will  be  applied  to  a double-body  problem. 

Let  (|)  be  the  nondimensional  perturbation  velocity 
potential  which  is  the  solution  of  the  double-body  problem. 
Then  (J)  satisfies  the  following  condition; 

4)_  = 0 on  z = 0 . (Al) 


Assume  the  Green  function  as  follows: 


G(P;Q)  r(P;Q)  r"(P;Q) 

where  r and  r'  have  the  same  definitions  as  in  Equation 


(31) . Then 


= 0 


z = 0 . 


Assume  the  same  coordinate  system  and  control  surfaces 
as  in  the  problem.  According  to  Green's  formula, 

4)  at  a point  P(x,y,z)  in  V can  be  written: 


())(P)  = 4:;jr  {()>(Q)  Gn(P;Q)  - <t)n  (Q)  G(P;Q)}dS.  (A3) 
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Using  the  conditions  (Al) , (A2)  and  following  the 

argument  as  in  Chapter  III,  the  equation  (A3)  can  be 
reduced  to ; 


<})(P) 


4it 


r 

. . 


{<|)(Q)  Gn 


H 


- (J)n(Q)  G}  dS 


+ 


47: 


fO 

dpr(p) 

-AR 


»00 

Gn  dC  . 

0 


(A4) 


Finally  <(>  on  H is  the  solution  of  an  integral  equation; 


^ 2tt 

<|>n  G)  dS 

H 

rO 

»00 

+ ^ 

Tir  J 

d pr(p) 

(A5) 

The  integral  equation  (A5)  is  identical  to  the  integral 
equation  for  problem  except  for  the  plus-minus  change 

in  the  Green  function.  Hence  the  mathematical  method  and 
the  computer  program  developed  for  solving  the  (ti^j  problem 
can  be  used  for  solving  the  double  body  problem  by  only 
changing  the  sign  in  the  Green  function.  The  solution  of 
the  double  body  problem  can  be  checked  by  various  experi- 
mental or  theoretical  results  of  wing  theory.  Hence  the 
theory,  the  computer  program  and,  indirectly,  the  solution 
of  the  problem  can  be  certified.  This  is  desirable 

since  the  (l)^  problem  is  rather  an  artificial  one  and 
hence  cannot  be  checked  by  comparing  with  experimental 
results . 

The  numerical  calculation  was  performed  for  the  strut 


shown  in  Figure  5.  A numerical  scheme  depending  on  the 
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finite  element  method  was  used.  In  Figure  A1  and  Figure 
A2  are  shown  the  distributions  of  nondimensional  pertur- 
bation velocity  potential  over  the  midsection  contour  for 
various  aspect  ratios.  In  symmetric  cases,  the  solutions 
should  coincide  with  the  two-dimensional  solution  on 
the  midsection  contour  for  any  aspect  ratio,  since  this  is 
a double-body  problem.  This  is  clearly  shown  in  Figure  Al. 

In  the  asymmetric  cases  shown  in  Figure  A2 , the  solutions 
calculated  approach  the  analytically  obtained  two-dimen- 
sional solution  as  the  aspect  ratio  increases.  The  solu- 
tions also  show  infinite  perturbation  velocity  along  the 
sharp  leading  edge.  The  calculated  side  force  coefficient 
is  plotted  vs.  aspect  ratio  in  Figure  A3.  The  side  force 
coefficients  for  rectangular  flat  plates  calculated  by 

j 

Daoud  (1973)  and  Kerwin  and  Openheim  (1974)  based  on  thin 

wing  theory  are  also  shown  in  Figure  A3  as  well  as  the  ; 

experimental  results  for  rectangular  wings  from  Thwaites 

i 

(1960:  p.  343) . 

From  the  results  shown  in  these  figures,  it  can  be 
concluded  that,  adopting  the  present  method,  reasonable  j, 

results  can  be  obtained  using  rather  a small  number  of 

i 

elements  in  a numerical  calculation.  It  seems  the  accuracy  i! 

i' 

of  the  present  numerical  solution  is  not  enough  for  dis- 
cussing the  effects  of  the  thickness  of  the  strut  or  the  ^ 

shape  of  the  bottom;  however,  the  accuracy  of  the  solution 


can  be  improved  to  any  required  degree  by  increasing  the 
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